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^ ■ Abstract 
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QQ , With an assumption that in the Yang-Mills Lagrangian, a left-handed fermion and 

m 

a right-handed fermion both expressed as the quaternion make an octonion which 
possesses the triality symmetry, I calculate the magnetic mass of the transverse self- 
dual gluon from three loop diagram, in which a heavy quark pair is created and two 



^ ■ self-dual gluons are interchanged. 

The magnetic mass of the transverse gluon depends on the mass of the pair created 
quarks, and in the case of charmed quark pair creation, the magnetic mass rrimag 
becomes approximately equal to Tc at T = Tc ~ 1.14A-jg^ ~ 260MeV. 

A possible time-like magnetic gluon mass from two self-dual gluon exchange is 
derived, and corrections in the B-meson weak decay vertices from the two self-dual 
gluon exchange are also evaluated. 
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1 Introduction 



In 1980, Linde[T] pointed out difficulties of infrared problems in the thermodynamics of 
mass-less Yang-Mills gas, and a possible solution via gluon acquiring the magnetic mass. At 
zero temperature, an effective gluon mass of 500 ± 200 MeV was predicted by Cornwall [2] . 
A problem that at relatively low temperature or high density, the pressure of the QCD 
thermodynamic gas becomes negative was pointed out in [3] and difficulties in the calculation 
of the thermodynamical potential in the order are discussed injH |5l [6]. 

The ground-state energy of the finite temperature quark-gluon system was classically 
derived from the Yang-Mills Lagrangian by Freedman and McLarren[7]. An extensive review 
of calculations done before 2004 are given in [8]. 

In QCD, the gluon is screened in the plasma through gluon loops, quark loops and ghost 

loops. In 1993 the non Abelian Debye screening mass was found to be sensitive to the 

non perturbative magnetic mass of the gluon[9J. By adding a magnetic mass vn^raag ^'^ 

transverse self-energy YlTiK), the instability could be evaded, when rrimag = (^f^d'^^cT and 

Cf > 1. rrimag IS cxpected to appear in the perturbative calculation in the order of and 

higher, but no systematic calculation of it was found. 

Alexanian and Nair[TU] calculated the magnetic mass of the gluon in a model based 

Cg'^T 

on Chern-Simons theory and obtained rrimag ~ (2.384) , where C = N for SU(N) 

47r 

gauge theory. There is, however, other models that yield different values of the magnetic 
massfTTj [12] . and the situation is worrisome. The best that one could do perturbatively to 
get the Debye mass isjl3j 




where N^, is the number of colors and Nj is the number of fermion flavors, but the value of 
c/ is left open. 

There is a systematic investigation of vacuum polarization tensor T{^y{kQ,k) of rela- 
tivistic plasma [14j. The collective plasma effects are characterized by the frequency Uk = 

2 



1 gT. The self-energy depend on four vector -u" of the fluid and the momentum i^" 
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of the virtual particle via two invariants uj = K'^Ua and k such that K°'Ka = cu^ — k"^. The 
transverse self-energy function YlTikjU) and the longitudinal self-energy function 11^(0;, /c) 
are calculated by a pair of four vector i?" and 5", each orthogonal to m", defined by 



When UJ > oOp, plane wave solution of -E" and i?" exist, but when u < Up they are screened. 
It was shown that the screening length for the magnetic fields diverges at w — ?■ 0, thus 5" is 
screened except at = 0. 

In [TB| [T7]. I proposed a calculation of the Domain Wall Fermion(DWF) propagator 
in quaternion basis and expressed vector fields in terms of the Pliicker coordinate, following 
E. Cartan[18]. In this framework, the spin structure of the quark pair in the self-energy 
diagram with two self-dual gluon exchange is uniquely defined when the color and spin of 
the incoming gluon are fixed. The process is similar to the one investigated in the technicolor 
theory [20]. 

In order to calculate 3°" at a; = 0, I adopt a model, in which a quark pair is created, and 
two self-dual gluons are exchanged and the pair is annihilated [TB| [T7] . From a left-handed 
quark described by a quaternion and a right-handed quark described by another quaternion, 
one can construct an octonion. An octonion posesses the triality symmetry [181 HH] and the 
spin structure of the magnetic coupling of a self-dual gluon to the quark loops can be fixed 
when the exchanged internal gluons are restricted to be self-dual. There is no direct evidence 
that the triality symmetry plays a role in the nature, but the difference of about a factor of 3 
in the effective flavor number for opening the conformal window in the Schrodinger functional 
scheme [21] v.s. in the momentum subtraction (MOM) scheme [22] could be understood, if the 
Schrodinger functional does not select a triality sector but MOM scheme selects one triality 
sector. Using the Banks-Zaks expansion |23]. Grunberg|25] showed that the non-perturbative 
effect modifies the critical flavor number of the perturbative QCD which was around 8[21j 
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to 4. Infrared stable QCD running coupling was observed also in holographic QCD[27], and 
in the polarized electron nucleon scattering [26 j. 

In order to check the relevance of self-dual two gluon exchange in the coupling ver- 
tex, I calculate the vertex of a B meson weak decay vertex into lepton and anti-neutrino. 
The decay propability of a B meson into a lepton and neutrino is measured by the Babar 
collaboration [31] and by the Belle collaboration |32]. Possible deviation from the standard 
model via Penguin diagram is discussed in [33] and I compare with the possible correction 
from two self-dual gluon exchange diagrams. 

The structure of this paper is as follows. In the sect. 2, I summarize the magentic mass 
problem and present the calculation of the self-energy using the quaternion bases. In sect. 3, 
the B-meson weak decay vertex including the two self-dual magnetic gluon exchange is 
investigated. Conclusion and discussion are given in the sect. 4. 

2 The self-energy of a gluon via exchange of self-dual 
gluons between heavy quark and heavy anti quark 

In 1993 the non abelian Debye screening mass was found to be sensitive to the non per- 
turbative magnetic mass of of the gluon[9j. The importance of the magnetic mass was 
recognized in [3], as he calculated the pressure of finite temperature plasma, which is de- 
fined a.s p = —TlogZ/V where Z is the partition function. The pressure derived at a 
high-temperature region 



becomes negative when extrapolated to low temperature T or to high chemical potential /i. 
The negative pressure implies that the convergence of the perturbation series breaks down, 
and a QCD phase transition at this point was discussed. 



P 
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Figure 1: The self-energy diagram of trans- Figure 2: The self-energy diagram of trans- 
verse polarized gluon through exchange of verse polarized gluon through exchange of 
self-dual gauge fields, flu self-dual gauge fields. II22 

In the infrared, however, the ring sum term log{l + IIt{K) / K'^) was found to contain 
infrared singular term, indepent of the choice of the gauge. In order to evade the equation 
of the one loop transverse gluon propagator near T = 

+ Uriko = 0, fc) = fc^ - g'^N.T^^^^^^^k = 

64 

be satisfied for positive k, a modification IIt{K) — )• IIt{K) + m^^^ was proposed [6J. With 
an arbitrary parameter c/ > 1, rrimag = Cf^g'^NcT ~ Cf^T was expected for as{k) = 0.3 
at /c ~ IGeV. The magnetic mass contributes to the pressure of the quark- gluon plasma in 
the order of g^T"^. In lattice simulations, the magnetic gluon propagator does not show a 
peak at the zero momentum and the pole structure of the magnetic mass was, excluded |35j. 
although how to detect the timelike pole on the lattice remains an open quest ion |34j. 

I do not assume a magnetic mass of g'^T for a producion of a g^ term in the pressure, 
but I consider a three loop diagram of order g^ including two self-dual gluon exchange, as a 
building block of the stabilized quark-gluon plasma system. The model does not contradict 
with the lattice simulation. 

In I discussed the quark propagator using quaternion [IS] basis. In this flamework, 

I consider tensor coupling of the external gluon field to a heavy quark internal self-dual vector 
field xi,X2, X3 couple to the fermion spinors. 

I consider the self-energy or the magnetic mass of a gluon, which is polarized in the 
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transverse direction. The transverse gluon field is 

where A" is the SU(3) color basis. The tensor coupling in the case of nm are, 74(7273 — 
7372)72 = 7i75 in the vertex and in the case of 11^22 is 74(7371 - 7i73)/2 = 7275 
The coupling of Fij{q) to the quark is 

y ■ 

= i'* {k)-f5eijn—(j'i^{p){eijnP'Aj{q) H ) 

= r{k)l,^^^^^{p){<y X A{q))n (2) 
m 

Internal self-dual gluon and the heavy quark couplings in and n22a are 

74(7374 - 7473)72 = -73, 74(7274 - 7472)72 = -72, and 74(7174 - 7173)72 = -7i- 
In IIii, I choose the quark represented by ^4 to be at rest and the self-dual gluon xz and 
X2 have momenta ky and p^, respectively. The quark ,^12 and ^31 have momenta —k and 
— p, respectively, and ^314 and ^124 have momenta Pz and ky, respectively. The propagator of 
quark ^23 have the numerator 73^^ -|- 72/^2/, but 73 is multiplied by 72 at the junction to the 
quark ,^314, and 72 is multiplied by 73 at the junction to the quark ^124 and effectively the 
numerator is proportional to as required by the assignment of ^23- 

I use the Clifford products rule (ao -l- a)(fco -\-b) — aobo + a^b -\- h^a — ah-\-axhto 
evaluate Chfford products of the following bases. 

C(t> = 6234 - 63* - 61J - C12fc 

Clp = 623 - 634* - ^314^' - il2Ak 

(t> = Co + ^14* + 64i + e34fe 

= iA + iii + i23 +iik (3) 

The xi gluon coupling to quark pair with 2 self-dual gluon exchange 11"^ in Coulomb 

gauge consists of two types i.e (— {12C314 + C31C124), or (^24^3 "^34^2)- The possible quark-anti 
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quark state between the self-dual gluon exchange of the former is ^4^23 and the latter is 

^1234^1- 

Since the trace of the two types are the same, I consider the amplitude 11°^ of the xi 
gluon 

Ja{k,p) = m g Tr7i75— — 5-737^ j72 — 5- ^ 

k^xpy'-f^k^ + m 'y2Py + l3kz + m i2Py + m 

X7175 7^- ^7275 — 2 ^1.2 A ^7375^- 2 ^ 

m k^ + Py + ki + Py + 

Similarly, the self-energy from 1122 becomes, by choosing the intermediate quark-anti 

quark state (-^3624 + ^12634), or (-^14^3 + ^346), 

Jb{k,p) = m g Tr-^2l^—— ^717^ ^73 — — 

m k^ + 171"^ k^m pj pj + 

k^xpz^ik^ + m izPz + likx + m jsPz + m 

X7275 ry— — ^7375 — 21/21 ^^i^s^^ ^ 

m k^ + m'^ Pi + k^ + p^ + m'^ 

The trace in the numerator of the Hn 

4(m^ + ^m\pl + kl) - m^plkl) (6) 
The numerator of the 1122 is 

A{m' + ?>m\pi + kl)-mYzkl) (7) 

I integrate numerically 

r-oo ^00 ^6 ^ 3^4(^2 ^ p2) _^2^2p2 ^^^^ 



lo Jo {k'^ + rn'^yip'^ + m'^y{k'^ + p'^ + m'^) {2'n'y 

where the factor 4 comes from integral on k and p from —00 to +00. 

I chosse the scale m = 1 and evaluate 

l + 3{x + y)-xy dxdy 
Jo Jo ix + my + inx + y + l)i27ry ^ ^ 

by varying A. Numerically, the integral is approximately 

At T = 0, when the wave function of the quark pair is available from lattice simulation, 

one can perform the integration over q numerically and obtain a non-screened magnetic mass. 
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At T 7^ 0, it is necessary to replace the self-dual gluon exchange propagator ——in 
J{k,p) by — 5 — 77^^- o — 77^ and check the self-consistency. 

In a quenched SU(2) lattice Landau gauge simulation, the temperature dependence of the 
transverse gluon and longitudinal gluon from temperature to twice the critical temperature 
Tc were recently measured |2 8 1 129] . Near the T = T^, the transverse gluon propagator showed 
a peak near q = 0.4GeV, but the longitudinal gluon propagator did not show this behavior. 
In a quenched SU(3) lattice Landau gauge simulation [30], the transverse gluon propagator, 
or the magnetic gluon propagator showed smooth momentum dependence in the range from 
T = 0.867Tc to T = 4.97Te. 

In the present work, I ignore these temperatures and the number of color dependence, 
and evaluate, using the bose Bose-Einstein distribution for the gluon, the magnetic mass at 
finite temperature T as 



m™,,(T)=T W-^tP^, (10) 



2 /(g) dQ 

Since q is in GeV, I use ks = 8.6 x 10~^eV/K= 8.6 x lO^^^GeV/K as the unit, and express 
the temperature in the unit of Tc = LMAj^^ [8J and Aj^s ~ 230MeV. 

There are two polarization directions and two different time orders, and assuming that 
the coupling of self-dual gluon and the external gluons to quarks are the same, I estimate the 
/(g) using the numerical results of T = as /(g) = g{qY — ^7 — 7^ = ces{q)^— — r— — . Here, 
I adopt as{q) obtained from Lattice simulation [22] which is consistent with the prediction 
of the holographic theory [27]. The parametrization of a{s) is [26] 

log I g^+"^g(g)^ 



A2 

where 

n{q) = vr(l + L , 2 , kJ, ^ /a^ + (^^)']"') (^2) 

log(m^/A^)(l + g/A) — 7 

and 7 = 1^ = 3!^, mg{q) = j^^, m = 1.024GeV, a = 3.008GeV-\ d = 0.840, b = 
1.425GeV-\ c = 0.908 and A = 0.349GeV. 
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Figure 3: The magnetic mass of the gluon through charm quark loops, as a function of T/Tc 
and a fit. 

The charmed quark-anti quark pair creating 3 loop diagram gives the magnetic mass of 
the order of near T = T^. The magnetic mass mmag{T)iii this case can be parameterized 
as 

mmagiT)/T = 0.711(T/T,)2 + 0.563(r/r,)" - 0.0628(T/T,)'^ 

The self-energy decreases as the quark mass of the loop increases. The magnetic mass 
of a gluon through bottom quark loops is less than fc^Tc and the thermal response of the 
bottom quark and the charmed quark differ qualitatively. 

When T = 0, I evaluate the self-energy as the component of the Clifford product. In the 
Coulomb gauge, the exchanged self-dual gluons are assumed to be xi,X2 and the quark-anti 
quark state that couple to can be ipCcp — )■ (^1.^23 + (,2C,3i), or (f)Cip — )■ (^14.^234 + ■^24^314)- 
The intermediate quark-anti quark state independent of the momentum could be ^1234^4 in 
the former case and ^123'Co in the latter case. However, since these states don't have virtual 
momentum excitations, I consider, instead transition of ipC^' to (pCip and vice versa in the 
intermediate state. 

I define self-energy U^^ for exchange of self-dual gluons Xi,X2 is obtained by using eigen 
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states (^14634 + 64^34) or (^23^1 + ^136)- 



J,, ^ i &^xrr -liK + m 1 -73(A;^ X g ) +m 
Jc[k,p) = m g — I r7 



m 



Py l2Py + m IsiPy X r^) + m li{qyPy)r^ + m 



7i75 



1 -'l2ir^k^)qy + m 



m 



(Py X r^)2 + m2 {qyPyYrl + 



(13) 





^234 ^12 ft 14 



Figure 4: The self-energy diagram of time-like gluon through exchange of self-dual gauge 
fields, 



1-44 



When the exchanged self-dual gluons are X1X2, the trace in the numerator of Jc{k,p) 
becomes 

4m^[(py^ + kj)qyrx - (1 + iqyrxy)kxPy][m'^ - {K x qy){Py x Tx)] (14) 

Numerically one could fix the momentum transfer {kxrx)qy —Py and integrate over Tx-, kx and 
Py, or the momentum transfer {qyPy)rx — kx and integrate over qy,Py and kx- 

When the external gluon is a spacelike photon, the infrared divergence can be regularized 
and yields anomalous magnetic moment | = 1 + ^- Since regularization of qy, Vx is necessary, 
the numerical calculation of the self-energy of the timelike gluons and photons remain as a 
future problem. 



3 The B-meson weak decay vertex 

When the two self-dual gluon exchange is important in the heavy quark heavy anti quark 
pair to a gluon coupling, similar two gluon exchange in B meson weak decay is expected to 
be important. I consider B — )• r^v decay, in which b quark u quark couple to W boson 
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through Cabibbo-Kobayashi-Maskawa matrix as 



{u{q + p), c{q + p), t(g + p))7m(1 - lb) 



\ 



Uud 


Uus 


Uub 






Ucb 


Utd 


Uts 


Utb 



( dig) ^ 



w^ipy + he 



= u{q + p)*747m(1 - lb)Uuh h{q) W^{p)^ + ■■■ 
and its hermitian conjugates. 



(15) 





Figure 5: The B — tv vertex. The direct 
term, xi is the transverse W boson. 



^ft 14^23 ^124 



Figure 6: The B — tv vertex. The three 
loop term. Xi is the transverse W boson, 
X2 and X3 are self-dual gluons. 



The direct coupling of W boson to the quark current has two types: 



Da{qy 



Tgw 
Tgw 



Tr- 



mt 



747i(l - 75) 



l2{qy +Py) + rriu 



-Akz{qy +py) 



dpydkz 



4>{Qy+Py,kz 



dk^dpy 
(27r)2 



{kl + ml){{qy + pyf + ml^ 



<P{,qy+Py,k^ 



(16) 



and 



Db{qz) = Tgw 
= Tgw 



-jsky + rrib 



Tr ' r" 747.(i - 75) T\ :::: ^i<i^ + 



k"^ + ml 
-4:ky{q^ + p^ 



{qz + PzY + ml 

dpzdky 



i2n) 



{kl + ml){{qz+Pzy + ml\ 



<P{qz+Pz,ky 



(2vr) 



(17) 



where (/)(p, k) is the B meson wave function in momentum space normalized as 

|0(p, k)\'^dpdk = 1. 



At large q, the Fourier transform of Gaussian wave function D{q) decreases too fast. 
Thus, I adopt the discrete cosine transform of the Bethe Salpeter equation evaluated by a 
lattice simulation |42]. Details of the choice of 0(p, k) is given in the Appendix. 
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The numerical value of Da{qy) at qy = ——fM^ ~ ~2.34 GeV is 0.6J\fTgw, where 
A/" = 0.025 is the wave function normahzation. At qy = 2.34GeV, the magnitude is about 
10~^ times smaller. 

In the 3 loop diagram, 1 choose the quark momentum is as the 6— quark. We are 
interested in the momentum region small compared to the temperature T, where the thermal 
part of the gluon propagator is simplified as [TO] 



- ^[^(^0 - u;,) + 6{ko + u,)] (18) 

where Uk = \/P~+~M^. 

The transition amplitude is obtained by calculating the trace, 

z 



s 73(^2 + kz) +muPyX -f2Py + Isikz + qz) +muPyX k^ 

X 1 - 75 7l74 . , L N2 I 2 7275 . ^2 , 2 I 2 ^3^5 

[qz + k-^Y + mi vfiu (fc^ + qzY + PZ + tHu 



l2Py + mu,, , .dpydk. 



l + ml i2nY 

. . _„2^2^Q^2^ _L ^3 _L _L ^ ^ _L l'^. _L ^ ^^2^ 



Ta 4 
gw9 



Aklpll'imlrriu + ml + k^{q^ + k^){-mb + rriu) + (mf, + m„)p2) 
{kl + ml){{q-, + k^Y + ml){pl + ml){pl + ml){{q^ + k,Y +Pl + ml) 



{kl + mljmpl + mL,(T))'^^^- (2vr)^ ^'"^ 



^ X 4 4 ffrj. -lik^ + rrib TIT -^sPz + 

X 7i(gx + A;^) +muPzX k^ 73^^ + 71 (g^ + k^) +m^p^x k^ 

x(l - 75)72747 , , X2 I ^ 7375 — 7T~yn — 2 ^i^s 

[qx + /i^x)^ + mu Pi + (g^ + k^Y + 

pl + ml ^' (27r)2 

Aklpl{?>mlmu + + A;^(g^ + A;^)(-mb + m„) + (m;, + mu)pl) 



Ti 4 



(A;2 + m2)((g^ + /c^)2 + m2)(p2 + m2)(p2 + m2)(p2 + (g^ + fcJ2 + ^2) 
+ (20) 



(fc^ + <a,(T))(pH<a,(T))"^"^' (27r)2 

When T = Tc, qz = —2.34 GeV, and the numerical value of the integral Ca{T,qz) is 
2.75A/'T^m„(7vi/(2Q;(ge//))^ < ■^T'^Qw- Here, g"^ is evaluated from phenomenological 27r(2Q;(gefr)), 
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which depends on the effective momentum transfer of self-dual gluons X2 and X3. I assumed 

In the case of B meson decay, the momentum transfer q is of the order of Mb/ 2, and the 
damping of (pip + q, k) etc are important, and for an exact evaluation, lattice simulations are 
necessary. 

There is no infrared divergence even when mmag{T) = 0, and the finite mmag{T) makes 
the vertex suppression at high T. 

The decay rate T = - oi B ^ W in PQCD is[36l EZ] 

T 

1 irp 2 

r(5 ^ = ^\Uu,\'GlftMl\M\\l - -j^f (21) 

where (0|'U7^756|i?^) = ifsPB = if^^(-'- ~ 1^)^- Including charged Higgs effect, 

fB\Uub\ ~ 0.85 is reported [36]. 

T{B ^ Iv) = -\U^,\'Ms\D{T) + C{Tni - (22) 

The correction C{T) to D{T) is not too large in the case of S„ = hu meson. Experi- 
mentally, t{B^) > t{B^) » t{Bc^). a Higgs enhancement effect in 6 rz/ + X [36] and 
R-parity violation in B ^ W [37] are discussed with an estimation /^^ = 450MeV[S7]. 

4 Discussion and Conclusion 

In this paper, an u independent, temperature dependent magnetic mass of a gluon is calcu- 
lated in a model with triality symmetry. I checked the consistency of the model by estimating 
a correction to the width of a B meson decay into a lepton and an anti-neutrino. The idea 
of triality symmetry comes from the difference of the critical flavor number for opening 
the conformal window in the Schrodinger functional method [2T] and in the Domain Wall 
Fermion method^!]. The polarized electron experiment at JLab[26] and the QCD running 
coupling constant in holographic QCD [23 and in the Coulomb gauge lattice simulation in 
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MOM scheme [22] suggest that the flavor number 2+1 system is not far from the conformal 

window, but the Schrodinger functional method suggests that the flavor number should be 
more than 8 for the conformal window opens. If DWF or MOM scheme selects one triality 
sector, while the Schrodinger functional does not due to the difference of the boundary con- 
dition of the wave function, the difference of about a factor of three in the critical number 
of flavors could be understood. 

The non-perturbative effects in finite temperature QCD is difficult. In [38], problems in 
renormalization of the three loop diagram as J{p,k) or C{p,k) are discussed. If I replace 
the propagator denoted by ^4 is replaced by a point, an integral like 



appears. The momenta of internally exchanged gluons p and k in my model correspond to 
their t and q. The momenta corresponding to ^23 or ,^31 in Figj2]or FigJTjin [38] are not fixed, 
since a coupling of the external gluon to the quark ^4 is considered in their model. Their p 
is taken as an integration variable, but this complication is absent in our model. 

I approximated the quark wave function by the Bethe Salpeter wave function derived 
in a lattice simulation [42]. It would be possible to refine the wave function and extend to 
the weak decay of B^, which is discussed in [37] as a check of the R-parity violation in the 
Minimal Supersymmetric Standard Model(MSSM). How the vertex correction of the present 
model affects the charged Higgs Boson effects in the weak decay of B meson [36] is left to the 
future. 

In the model of domain wall fermion, the quark spinor possesses an arbitrary f/(l) phase 
that one can adjust to one triality sector of the lepton in the detector. I expect that there 
is three fold U{1) symmetry in the quark sector, but the symmetry is broken in the lepton 
sector and one triality is selected. The situation seen from the side of a quark is similar to 
the color-flavor locking[39]. In color-flavor locked states, the photon coupling to quark flavor 
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7 and gluon coupling to quark color F have the same structure 



2/3 







7/r 







1/3 



y -1/3 j 



and the combination 7 



^7 defines a massless gauge boson. 



This hypothesis might explain the origin of dark matter or the weakly interacting massive 
particles (WIMPs), which are not detected by electro- magnetic detectors on the earth. If 
the vector field of a photon can also be expressed by using Pliicker coordinates, it could also 
acquire a magnetic mass. 
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Appendix: Quark wave function (f){p, k) 



The discrete cosine transform(DCT) is used in image compressions etc. For a given data 

Tik{2x + 1) 

s{x) (x = 0, 1, ■ ■ ■ — 1) one multiplies cos and make a weighted summation 

. When the DCT is applied to \ —e~^^'^^'^ , an expresion close to its Fourier transform 



TT 



' — -e ^^'^^^ is obtained. 
I consider the case (3 = 2 and evaluate DCT of 4.9e~^ . The factor 4.9 is the normalization 
that is used in the Bethe Salpeter wave function of the B meson which will be discussed later. 
When the DCT is applied, the scale of ordinate changes by the factor and that of abscissa 
changes by 27r. The DCT scaled by multipying 1/2 is compared with the original in FigJTl 
When the maximum of a; = 2 corresponds to Ifm, Ax = 0.05 and in momentum space 
Ap = 2n fm-i=1.18 GeV. 

The parameter of the Gaussian wave function is taken from the flux-tube model [IHl 
1^ . The model was applied to low energy meson decays and meson-baryon couplings and 
produced reasonable results. 



Figure 7: The discretized gaussian wave Figure 8: The discretized Bethe Salpeter 
function and its DCT. wave function. 



The Bethe Salpeter wave function of the Qq meson state was calculated in [32]. They 
adopted the scale ro = 3a = 2.68GeV^^=0.529fm, and the s-wave w{R) can be parameterized 



as 



iu{ — ] = 4.9e '°'^^-o = 4.9 X 7 'o , 
To, 
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The wave lunction is normalized as / w \ — — — ~ 1. 

/o VoJ VoJ To 

R 



I take 20 points of ^ from to 1.9, which corresponds to radial distance between h anti 



quark and u quark of about Ifm. I make a discrete cosine transform of this wave function 
shown in FigJHl I multiply the scale of the ordinate 1/2 at a moment. Since AP corresponds 
to 27r/(2 X 0.5287)fm~^ = 5.94fm~^ = l.lTGeV, I change the scale of abscissa and plot in 
Figjni The wave function w{k) obtained by DCT is normalized as 

19 . 



0=0 



10 



When the scale of the ordinate is multiplied by \ as in the case of Gaussian, and abscissa is 
transformed to GeV, 

19 

^*(jr47r(1.17(j + l)GeYf x AP ~ 



^-[(J^)2_(J^)2]/(2/3) 

and the Bethe Salpeter wave function corresponds to the e part, I identify —w{j) 



Since the Gaussian wave function can be decomposed as e 

uhe Bet 
fP — 



as 



V2 

infinitely heavy. 



and the wave function of the center of mass to be 1, i.e. I regard h quark is 




Figure 9: The DCT of the discretized Bethe Salpeter wave function and the DCT of the 
Gaussian wave function. 



The wavefunction k) is approximated as a function of \p — k\/ \/2. 
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p 





2tt 


47r 


Gvr 


87r 


fm ^ 


4^ gauss 


5.10 


3.80 


0.44 


-0.01 


-0.04 




P 





5.94 


11.9 


17.8 


23.8 




4>MP 


3.03 


2.70 


1.19 


0.65 


0.37 





Table 1: The momentum dependence of the relative wavefunction of quarks. 
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